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The  main  focus  of  this  paper  is  a  study  that  empowers  us  to  understand  how 
the  temperature  variation  affects  the  transmission  line  resistance  and  as  a 
result  the  power  flow  analysis  with  a  specific  end  goal  to  assess  losses  in  the 
electrical  network.  The  paper  is  composed  of  two  sections;  the  first  part  is  a 
power  flow  study  under  normal  conditions  utilizing  the  neural  network 
approach  while  the  second  investigated  extreme  learning  machine  algorithm 
efficiency  and  exactitude.  Extreme  learning  machine  algorithm  has  been 
used  to  settle  several  complications  in  power  system:  load  forecasting,  fault 
diagnosis,  economic  dispatch,  security,  transient  stability;  Thus,  we  proposed 
to  study  this  technique  to  figure  out  this  sort  of  complex  issue.  The  study 
was  conducted  for  IEEE  30  bus  test  system.  The  simulation  results  are 
exposed  and  analyzed  in  detail  at  the  end  of  this  paper. 
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1.  INTRODUCTION 

Eleetrie  power  systems  are  real-time  energy  delivery  systems.  Real-time  means  that  power  is 
generated  transported  and  supplied  the  moment  we  turn  on  the  light  switeh  [1].  Eleetrie  power  systems  are 
not  storage  systems  like  water  and  gas  systems.  Instead,  generators  produee  the  energy  as  the  demand 
ealls  for  it  [2,  3]. 

The  huge  amount  of  power  generated  in  a  power  station  (hundreds  of  MW)  is  to  be  transported  over 
a  long  distanee  (hundreds  of  kilometers)  to  load  eenters  to  eater  power  to  eonsumers  with  the  help  of  high- 
voltage  transmission  line  in  order  to  reduee  transportation  losses.  Transmission  and  distribution  losses  reaeh 
a  rage  of  4%  to  5%  of  the  total  eonsumption.  The  higher- voltage  transmission  lines  require  larger  struetures 
with  longer  insulator  strings  in  order  to  have  greater  air  gaps  and  needed  insulation. 

Conduetor  material  (all  wires),  type,  size,  and  eurrent  rating  are  key  faetors  in  determining  the 
power  handling  eapability  of  transmission  lines.  A  eonduetor  heats  up  when  eurrent  flows  through  it  due  to 
its  resistanee,  eonduetors  tend  to  inerease  their  resistanee  with  an  inerease  in  temperature,  rated  by  how  mueh 
eurrent  eauses  them  to  heat  up  to  a  predetermined  amount  of  degrees  above  ambient  temperature.  The 
amount  of  temperature  rise  above  ambient  (i.e.,  when  no  eurrent  flows)  determines  the  eurrent  rating  of  a 
eonduetor.  For  example,  when  a  eonduetor  reaehes  70°C  above  ambient,  the  eonduetor  is  said  to  be  at  full 
load  rating  [1]. 

Seientists  used  to  mention  the  dynamie  thermal  rating  in  literature.  Conduetors  transporting  energy 
in  low  temperatures  are  more  effieient.  For  an  overhead  transmission  line,  the  maximum  temperature  it 
supports  is  80C°  [4-6].  The  eonduetor  transmitting  700  amperes  at  40C;  ean  transport  1200amperes  at  10C°, 
whieh  implies,  from  summer  season  to  winter;  transmission  line  eapaeity  ean  reaeh  the  double  [4,  6,  7]. 
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The  conductor  temperature  varies  in  the  range  of  0°  to  60°C  or  even  larger,  resulting  changes  of  the  line 
resistance  [3,  8,  9]. 

Meanwhile,  the  resistance  fluctuations  are  explicitly  taken  into  account  and  coupled  to  power  flow 
analysis.  Power  flow  fluctuations  can  cause  various  phenomena  ranging  from  frequency  and  voltage  stability. 
Thermal  contingency,  to  critical  scenarios  such  as  components  outages  and  black  out  [10].  For  the  purposes 
of  load  flow  determination  and  system  state  estimation,  a  variety  of  methods  has  been  used.  The  development 
of  these  methods  is  mainly  led  by  the  basic  requirements  of  load  flow  calculation,  which 
can  be  summed  up  as 

a.  the  convergence  properties 

b.  the  computing  efficiency  and  memory  requirements 

c.  the  convenience  and  flexibility  of  the  implementation 

Mathematically,  the  load  flow  problem  is  a  problem  of  solving  a  system  of  nonlinear  algebraic 
equations.  With  the  scale  of  power  system  continually  expanding,  the  dimension  of  load  flow  equations  now 
becomes  very  high,  thus  research  is  still  very  active.  Newton  Raphson  technique  stands  to  be  the  most 
important  approach  to  solve  the  load  flow.  Moreover  this  accurate  technique  became  unacceptable  for  real¬ 
time  electric  power  system  operation  (dynamic),  the  problems  result  mainly  from  the  speed  of  operation  and 
the  dependability  (no  missing  operations)  ,hence  the  introduction  of  the  artificial  intelligence  or  more 
precising  extreme  learning  machine  due  to  its  extremely  fast  training,  good  generalization,  and  universal 
approximation  /classification  capability  [11].  It  has  been  demonstrated  to  have  excellent  learning  accuracy 
/speed  in  various  applications,  such  as  face  classification,  image  segmentation,  human  action  recognition, 
prediction,  and  functions  approximation.  ELM  offers  effective  solutions  for  the  electric  power  system  as  it  is 
demonstrated;  in  the  article  [12]  a  novel  version  of  ELM,  which  can  be  applied  to  train  a  single  hidden  layer 
recurrent  neural  network.  ELM  was  presented  and  employed  in  order  to  achieve  higher  forecasting  electricity 
load  accuracy.  In  [13]  authors  introduced  an  ELM-based  predictor,  then  it  was  developed  for  the  power 
system.  The  predictor  has  been  used  for  real-time  to  enhance  the  dynamic  security  of  the  power  systems. 
Article  [14]  explores  a  fault  diagnosis  method  for  the  key  mechanical  components  of  the  wind  turbine.  The 
method  based  on  improved  extreme  learning  machine  (lELM)  combined  the  fault  characteristics  of  gear  and 
bearing  of  wind  turbine  drive  system.  A  robust  observer  to  detect  and  isolate  sensors  fault  in  a  variable  speed 
wind  turbine  is  discussed  in  [15].  The  designed  approach  was  based  on  the  ELM  algorithm  for  the  prediction 
of  unknown  model  dynamics  to  design  the  robust  non-linear  observer.  The  article  [16]  envisages  the 
application  of  craziness-based  PSO  to  enhance  the  performance  of  ELM  for  feature  classification  of  single 
and  combined  power  quality  disturbances.  Thus  we  introduced  this  fast  and  accurate  approach  in  this  paper  to 
examine  the  variation  of  the  transported  power,  load  flow  and  losses  due  to  the  influence  of  temperature  on 
conductor  resistance.  To  start  with,  we  attempted  to  study  the  first  proposition  which  is  the  neural  network 
Then  for  further  optimization  and  in  terms  of  precision,  exactness,  and  proficiency  extreme  learning  machine 
algorithm  was  applied  to  address  the  above  problem. 


2.  RESEARCH  METHOD 

2.1.  Neural  networks 

The  origin  of  artificial  neural  networks  comes  from  the  biological  neuron  modeling  test  by  Warren 
McCulluch  and  Walter  Pitts.  Initially  Neural  networks  objective  was:  patterns  recognition,  classification  then 
it  becomes  very  interesting  in  all  domains  [17]. 

2.1.1.  Neuron  Model 

Indeed;  the  simple  architecture  of  a  single  neuron  with  a  single  layer  is  composed  essentially: 

First,  the  scalar  input  p  which  is  multiplied  by  the  scalar  weight  w  to  form  the  product  wp  this  operation  takes 
place  at  the  integrator  which  is  one  of  the  important  parts  of  the  neuron;  the  result  is  again  a  scalar  added  to  a 
scalar  bias  b  to  form  the  net  input  n. 

The  result  n  is  then  transformed  by  a  transfer  function  /  which  produces  the  neuron  output  a.  [18, 
19],  the  model  is  shown  in  Figure  1. 
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Figure  1.  Artificial  neuron  model 


n  =  1.%!  w^jp  -  b 


(1) 


2.1.2.  Neural  Network  Learning 

Learning  is  a  process  by  which  the  free  parameters  from  a  neural  network  are  adapted,  through  a 
simulation  process  by  the  environment  in  which  the  network  is  contained  [17,  19,  20] 

The  supervised  and  unsupervised  methods  are  essentially  the  two  types  of  learning  for  a  neural  network.  In 
our  case,  we  are  going  to  use  the  supervised  learning  where  the  training  is  controlled  by  an  external  agent 
who  watches  the  answer  that  the  network  supposed  to  generate  from  the  determined  entrance.  The  supervisor 
compares  the  output  of  the  network  with  the  expected  one  and  determines  the  amount  of  modifications  to  be 
made  on  the  weights  [18,  21,  22]. 

The  supervised  learning  can  be  done  through  the  three  following  Paradigms: 

a.  Error  correction  learning 

b.  Reinforcement  learning 

c.  Stochastic  learning. 

During  the  artificial  neural  network  training  by  error  correction,  the  weights  of  the  communication 
links  are  adjusted  trying  to  minimize  a  function  of  cost  depending  on  the  difference  between  the  desired 
values  and  the  obtained  one  from  the  output  network  [19,  23]. 

It  is  required  to  determine  the  modifications  on  weights  using  the  committed  error,  the  correction  error  is 
calculated  as  the  following  expression: 

=  j?yy(4-yk)  (2) 

Where: 

dWfci is  the  weight  variation  of  the  connection  between  the  neurons  j  from  the  anterior  layer  and  the  output 
layer  node  k. 

y^\  the  desired  neuron  output. 

yy,  and  y/^  are  the  output  values  produced  in  the  neuron  i  and  k  respectively. 

7]  a  positive  factor  denominated  the  learning  rateO  <  77  <  1  [18,  24,  25]. 

We  can  say  that  the  new  weight  is  proportional  to  the  committed  error  produced  by  the  network  and  the 
answer  of  the  node  j  in  the  previous  layer  as  shown  in  the  following  formula: 

^actual  ^  ^anterior  ^3^ 


The  total  error  is  known  as  the  mean  squared  error  expressed  in  the  equation  (4)  [26,  27]: 

Error  =  —'Zp=i'Zk=i  {ylf^  ~  (4) 

Where: 

N  is  the  number  of  neurons  in  the  output  layer  and  P  is  the  number  of  examples  in  the  training  sample. 
Normally  to  minimize  this  criterion  function,  is  employed  learning  rule  given  by  the  descending  gradient  [18, 
19]. 


AWfcy  =  -  r?Vw|,j  (Errorgiobai) 


(5) 
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If  (  T])  takes  a  little  value  the  learning  process  is  made  smoothly  which  gives;  as  a  result,  an  increment  in  the 
time  of  convergence  to  a  stable  solution.  In  the  other  hand  if  ryhas  an  important  value  the  speed  of  learning 
gets  increased  but  there  will  be  a  divergence  risk  and  this  causes  the  instability  of  the  system  [25-27]. 


2.2.  Extreme  Learning  Machine  ELM 

Huang  et  al.  proposed  an  extreme  learning  machine  algorithm,  as  a  new  learning  scheme  working 
for  single  hidden  layer  feed-forward  networks  (SLFNs)  [28,29]. ELM  is  entirely  different  from  traditional 
iterative  learning  algorithms  as  it  randomly  selects  the  input  weights  (W)  and  biases  (b)  for  hidden  nodes  and 
analytically  calculates  the  output  weights  (p)  by  finding  the  least  square  solution. 

In  doing  so,  it  is  proven  that  the  training  error  can  still  be  minimized  with  even  better  generalization 
performance  [30].  Notably,  the  learning  speed  of  ELM  can  be  thousands  of  times  faster  than  traditional  feed¬ 
forward  network  learning  algorithms  like  BP  algorithm.  Compared  with  conventional  learning  algorithms, 
ELM  not  only  tends  to  reach  the  smallest  training  error  but  also  to  obtain  the  smallest  norm  of  weights  [29, 
31]. 

For  N  arbitrary  distinct  samples  (Xj,  tj )  where: 

(^il  >  ^i2  > . >  ^in  )  ^  R  (6) 

=  . ttmV  (7) 


Standard  SFLNs  with  Af~nodes  and  activation  function  g(x)  are  mathematically  modeled  as: 


lf=iPi3iWi.Xj  +  bi)  =  Oj 


(8) 


Where  Wj  =  (w^,  Wi2, . ,  ^  is  the  weight  vector  connection  (  ith)  hidden  node  and  input  nodes. 

I^i  —  . the  weight  vector  connecting  the  (ith)  hidden  node  and  output  nodes,  and  bf.  is 

the  threshold  of  the  (ith)  hidden  node.  Wj.Xy  denotes  the  inner  product  of  Wi  and  Xy. 

That  standard  SFLN  withAf~  hidden  nodes  and  ^(x)  activation  function  can  approximate  these  N  samples 
with  zero. 

i:p|o,-t,||  =  o 


What  means  there  exist  Pi,  Wi  and  bi  such  that: 

i:f=i  A5i(^y)  =  YJi=iPig(wi.Xj  +  bi)  =  tj 


Equation  (10)  can  be  written  compactly  as 
Hp  =  T 

Where  : 

//(wi  ...w^,Xi  ...x^,  hi  = 


g(wi.xt  +bt) 
.5(Wi.Xjv  +fti) 


5(w/v“-^i  +  V) 


r/^ii 

T 

p  = 

*** 

andr  = 

N~*m 

Xn- 

(9) 

(10) 

(11) 

(12) 


Thus,  the  extreme  learning  machine  algorithm  can  be  summarized  as  follows: 

Step  1 :  Randomly  assign  input  weight  VKiand  bias  bi,  (i  =  1,2,.,N). 

Step2:  Calculate  the  hidden  layer  output  matrix  H. 

Step3:  Calculate  the  output  weight  /?  by  P  =  H’*'T 

If  hidden  layer  output  matrix  is  non-square  (that  means  <  N),  the  determinant  of  H  det(H)  =  0  and 
(rank  H<AF^),  in  linear  algebra  matrix,  is  assessed  as  a  singular  matrix  that  cannot  be  directly  inverted  [32- 
34]. 


p=(H^H)-iH^T  (13) 

The  inverse  of  H  cannot  be  determined  if  H  is  not  a  full-rank  matrix.  In  this  case,  the  generalized 
inverse  matrix,  known  as  pseudo  inverse,  is  used.  In  order  to  calculate  the  pseudo  inverse,  various  approaches 


Power  flow  variation  based  on  extreme  learning  machine  algorithm  in  power  system  (Labed  Imen) 


1248  □ 


ISSN:  2088-8694 


can  be  used.  The  most  widely  known  type  is  the  Moore-Penrose  pseudo-inverse  [35].  The  pseudo-inverse  of 
the  matrix  H,  namely  can  be  computed  via  least  square  solution: 

//t  =  (14) 

This  solution  is  aeeurate  as  long  as  the  square  matrix  is  invertible.  However,  it  is  singular  in 

some  applieations.  In  SLFN,  it  is  singular  in  most  of  the  eases  sinee  there  is  a  tendeney  to  seleet  m  <<  n.  In 
ELM,  Huang  et  al.,  have  solved  this  problem  using  Moore-  Penrose  pseudo  inverse  of  H  matrix  to  overeome 
this  problem  [32,  36,  37]. 


3.  RESULTS  AND  ANALYSIS 

In  order  to  verify  the  proposed  two  approaehes;  ELM  and  neural  networks,  the  paper  study  is 
eondueted  on  IEEE  30  bus  test  system  shown  in  Figure  2,  it  eomprises  of  6  generators,  bus  number  one  (1)  is 
eonsidered  as  a  slaek  bus.  Likewise,  24  loads  powered  over  41  transmission  lines. 


Figure  2.  Single  Line  Diagram  of  IEEE-3 Obus  [38] 


3.1.  Line  Flow  Solution  Using  Newton-Raphson  Algorithm 

Table  1  displayed  the  line  flow  results  at  25  C°.  It  is  the  initial  step  to  provide,  aetive  and  reaetive 
power  and  the  total  system  lossesThe  whole  system  is  operating  under  normal  eonditions.  These  results  are 
used  to  be  the  input  matrix  for  the  neural  network  training  phase. 


Table  1.  The  line  flow  and  losses  at  T=25C° 


i  to  j 

P 

(MW) 

Q 

(MVAR) 

S 

(MVA) 

Losses 

(MW) 

Losses 

(Mvar) 

itoj 

P 

(MW) 

Q 

(MVAR) 

S 

(MVA) 

Losses 

(MW) 

Losses 

(Mvar) 

1-2 

177,778 

-22,148 

179,152 

5,464 

10,524 

12-  13 

0,000 

-10,291 

10,291 

0,000 

0,133 

1-3 

83,221 

5,127 

83,378 

2,808 

7,085 

12-14 

7,856 

2,442 

8,227 

0,075 

0,155 

2-4 

45,712 

2,705 

45,792 

1,106 

-0,517 

12-15 

17,857 

6,947 

19,161 

0,217 

0,428 

2-5 

82,990 

1,703 

83,008 

2,995 

8,178 

12-16 

7,208 

3,363 

7,954 

0,053 

0,112 

2-6 

61,912 

-0,958 

61,920 

2,048 

2,264 

14-  15 

1,582 

0,687 

1,724 

0,006 

0,005 

3-4 

78,012 

-3,158 

78,076 

0,771 

1,344 

15-  18 

6,014 

1,744 

6,262 

0,039 

0,080 

4-6 

70,126 

-17,526 

72,282 

0,604 

1,179 

15-23 

5,001 

2,956 

5,810 

0,031 

0,063 

4-12 

44,121 

14,646 

46,489 

0,000 

4,685 

16-17 

3,654 

1,451 

3,932 

0,012 

0,027 
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itoj 

P 

(MW) 

Q 

(MVAR) 

S 

(MVA) 

Losses 

(MW) 

Losses 

(Mvar) 

itoj 

P 

(MW) 

Q 

(MVAR) 

S 

(MVA) 

Losses 

(MW) 

Losses 

(Mvar) 

5-7 

-14,205 

10,500 

17,664 

0,151 

-1,687 

18-  19 

2,775 

0,765 

2,879 

0,005 

0,010 

6-7 

37,523 

-1,885 

37,570 

0,367 

-0,598 

19-20 

-6,730 

-2,645 

7,231 

0,017 

0,034 

6-8 

29,528 

-3,754 

29,766 

0,103 

-0,558 

21-22 

-1,877 

-1,594 

2,462 

0,001 

0,001 

6-9 

27,693 

-7,322 

28,644 

0,000 

1,594 

22-  24 

5,654 

2,788 

6,304 

0,043 

0,067 

6-10 

15,823 

0,653 

15,836 

0,000 

1,278 

23-24 

1,770 

1,293 

2,192 

0,006 

0,012 

6-28 

18,819 

-9,618 

21,134 

0,060 

-13,086 

24  -25 

-1,325 

1,602 

2,079 

0,008 

0,014 

8-  28 

-0,575 

-2,370 

2,438 

0,000 

-4,368 

25  -26 

3,545 

2,366 

4,262 

0,045 

0,066 

9-  11 

0,000 

-15,657 

15,657 

0,000 

0,462 

25-27 

-4,877 

-0,778 

4,939 

0,026 

0,049 

9-  10 

27,693 

6,741 

28,501 

0,000 

0,809 

27-28 

-18,184 

-4,157 

18,653 

-0,000 

1,309 

10-20 

9,027 

3,560 

9,704 

0,081 

0,180 

27  -29 

6,189 

1,668 

6,410 

0,086 

0,162 

10-  17 

5,372 

4,414 

6,953 

0,014 

0,037 

27-  30 

7,091 

1,661 

7,283 

0,161 

0,304 

10-21 

15,733 

9,842 

18,558 

0,110 

0,237 

29-  30 

3,704 

0,606 

3,753 

0,033 

0,063 

10-22 

7,583 

4,490 

8,813 

0,052 

0,107 

Total  loss 

17,599 

22,244 

3.2.  The  Neural  Network 

3.2.1.  The  training  phase  using  Levenberg-  Marquard  algorithm 

a.  The  training  algorithm  is  Levenberg-Marquard  (trainlm),  as  it  is  shown  in  (Figure  3). 

b.  The  performance  is  the  mean  squared  error  (mse)  (Figure  4). 

c.  The  convergence  is  satisfied  at  t=15sec  (Figure  3). 

d.  The  performance  value  reached  3. 1 5 e- 19  (Figure  4). 

In  tables  1  and  2  and  also  (Figure  5),  it  is  evident  that  the  neural  network  outputs  and  targets  are 
completely  identical.  The  main  squared  error  reaches  the  value  3.15e-19  (Figure  4).  Bias  and  weights  are 
adjusted. 

Here  below  we  can  see  clearly  the  neural  network  architecture  and  training  phase  exposed  in  Figures  3,  the 
performance  in  Figure  4  and  the  regression  in  Figure  5  respectively. 


Figure  3.  N.N  architecture  and 
training  phase 


Figure  4.  Performance 


Figure  5.  Regression 


Table2.  Fine  flow  based  on  a  neural  network  at  25C° 


itoj 

P 

(MW) 

Q 

(MVAR) 

Losses 

(MW) 

Losses 

(Mvar) 

itoj 

P 

(MW) 

Q 

(MVAR) 

Losses 

(MW) 

Losses 

(Mvar) 

1-2 

177,778 

-22,148 

5,464 

10,524 

12-  13 

1,705  e-13 

-10,291 

l,421e-14 

0,133 

1-3 

83,221 

5,127 

2,808 

7,085 

12-14 

7,856 

2,442 

0,075 

0,155 

2-4 

45,712 

2,705 

1,106 

-0,517 

12-15 

17,857 

6,947 

0,217 

0,428 

2-5 

82,990 

1,703 

2,995 

8,178 

12-16 

7,208 

3,363 

0,053 

0,112 

2-6 

61,912 

-0,958 

2,048 

2,264 

14-  15 

1,582 

0,687 

0,006 

0,005 

3-4 

78,012 

-3,158 

0,771 

1,344 

15-  18 

6,014 

1,744 

0,039 

0,080 

4-6 

70,126 

-17,526 

0,604 

1,179 

15-  23 

5,001 

2,956 

0,031 

0,063 

4-  12 

44,121 

14,646 

7,105e-15 

4,685 

16-17 

3,654 

1,451 

0,012 

0,027 

5-7 

-14,205 

10,500 

0,151 

-1,687 

18-  19 

2,775 

0,765 

0,005 

0,010 

6-7 

37,523 

-1,885 

0,367 

-0,598 

19-20 

-6,730 

-2,645 

0,017 

0,034 

6-8 

29,528 

-3,754 

0,103 

-0,558 

21-22 

-1,877 

-1,594 

0,001 

0,001 

6-9 

27,693 

-7,322 

-7,105e-14 

1,594 

22-  24 

5,654 

2,788 

0,043 

0,067 

6-10 

15,823 

0,653 

-7.993e-15 

1,278 

23-24 

1,770 

1,293 

0,006 

0,012 
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i  to  j 

P 

(MW) 

Q 

(MVAR) 

Losses 

(MW) 

Losses 

(Mvar) 

itoj 

P 

(MW) 

Q 

(MVAR) 

Losses 

(MW) 

Losses 

(Mvar) 

6-28 

18,819 

-9,618 

0,060 

-13,086 

24  -25 

-1,325 

1,602 

0,008 

0,014 

8-  28 

-0,575 

-2,370 

8,881e-14 

-4,368 

25  -26 

3,545 

2,366 

0,045 

0,066 

9-  11 

-2.273e-13 

-15,657 

-1.421e-14 

0,462 

25-  27 

-4,877 

-0,778 

0,026 

0,049 

9-  10 

27,693 

6,741 

0,000 

0,809 

27-  28 

-18,184 

-4,157 

-l,421e-14 

1,309 

10-20 

9,027 

3,560 

0,081 

0,180 

27  -29 

6,189 

1,668 

0,086 

0,162 

10-  17 

5,372 

4,414 

0,014 

0,037 

27-  30 

7,091 

1,661 

0,161 

0,304 

10-21 

15,733 

9,842 

0,110 

0,237 

29-  30 

3,704 

0,606 

0,033 

0,063 

10-22 

7,583 

4,490 

0,052 

0,107 

Total  loss 

17,599 

22,244 

The  NN  is  trained  and  ready  to  be  utilized  in  real  time  for  resistanee  variation  Rt=3oc°  and  Rt=o5c°. 
Line  flow  results  using  the  neural  network  approaeh  are  exposed  in  Table  3  for  the  two  resistanee  Values 

Rt=3oc°  and  Rt=o5C°. 


Table  3.  Line  flow  using  the  trained  neural  network  in  real  time  (T=30C°  and  T=5C°  ) 


itoj 

R(pu) 

30C° 

P 

(MW) 

Q 

(MVar) 

Loss 

(MW) 

Loss 

(Mvar) 

R(pu) 

05C° 

P 

(MW) 

Q 

(MVar) 

Loss 

(MW) 

Loss 

(Mvar) 

1-2 

0.0195 

178,010 

-23,312 

5,593 

10,591 

0,0177 

176,840 

-17,501 

4,960 

10,273 

1-3 

0,0460 

83,370 

4,788 

2,871 

7,119 

0,0416 

82,639 

6,479 

2,561 

6,961 

2-4 

0,0581 

45,709 

2,525 

1,127 

-0,520 

0,0525 

45,719 

3,432 

1,024 

-0,507 

2-5 

0,0481 

83,105 

1,320 

3,061 

8,208 

0,0435 

82,543 

3,232 

2,738 

8,069 

2-6 

0,0592 

61,909 

-1,262 

2,087 

2,263 

0,0535 

61,918 

0,268 

1,890 

2,270 

3-4 

0,0134 

78,099 

-3,531 

0,788 

1,350 

0,0121 

77,678 

-1,682 

0,703 

1,320 

4-6 

0,0121 

70,138 

-18,018 

0,618 

1,188 

0,0109 

70,080 

-15,564 

0,549 

1,147 

4-  12 

0,0000 

44,155 

14,581 

-  7,5836e-ll 

4,689 

0,000 

43,989 

14,901 

-l,660e-ll 

4,669 

5-7 

0,0468 

-14,155 

10,720 

0,156 

-1,682 

0,0424 

-14,395 

9,613 

0,133 

-1,704 

6-7 

0,0272 

37,485 

-2,102 

0,374 

-0,600 

0,0246 

37,669 

-1,010 

0,341 

-0,594 

6-8 

0,0122 

29,531 

-4,254 

0,106 

-0,556 

0,0110 

29,520 

-1,763 

0,094 

-0,563 

6-9 

0,0000 

27,693 

-7,341 

4,263e-13 

1,595 

0,0000 

27,690 

-7,249 

-l,070e-ll 

1,590 

6-10 

0,0000 

15,822 

0,651 

-3,552e-14 

1,2790 

0,0000 

15,823 

0,658 

5,770e-12 

1,277 

6-28 

0,0172 

18,809 

-9,686 

0,061 

-13,082 

0,0155 

18,857 

-9,348 

0,055 

-13,102 

8-28 

0,0648 

-0,575 

-2,286 

-7,816e-14 

-4,3680 

0,0586 

-0,574 

-2,702 

2,680e-12 

-4,371 

9-  10 

0,0000 

27,693 

6,772 

-5,115e-13 

0,8100 

0,0000 

27,690 

6,617 

6,050e-12 

0,8060 

9-  11 

0,0000 

l,421e-12 

-15,708 

-7,816e-14 

0,4650 

0,0000 

7,069e-13 

-15,456 

-l,584e-12 

0,449 

10-17 

0,0330 

5,361 

4,431 

0,015 

0,0370 

0,0862 

5,415 

4,347 

0,013 

0,037 

10-  20 

0,0954 

9,0270 

3,5720 

0,0820 

0,1810 

0,0298 

9,029 

3,509 

0,074 

0,180 

10-21 

0,0354 

15,7410 

9,8430 

0,1120 

0,2370 

0,0320 

15,701 

9,837 

0,101 

0,236 

10-22 

0,0741 

7,5870 

4,4880 

0,0530 

0,1070 

0,0670 

7,568 

4,499 

0,048 

0,106 

12-  13 

0,0000 

3,4104e-ll 

-10,3200 

l,6667e-09 

0,1330 

0,0000 

3,471e-12 

-10,174 

3,064e-10 

0,130 

12-14 

0,1255 

7,858 

2,433 

0,0760 

0,1550 

0,1135 

7,849 

2,474 

0,069 

0,155 

12-15 

0,0674 

17,8760 

6,9330 

0,222 

0,4290 

0,0610 

17,783 

7,002 

0,199 

0,426 

12-16 

0,0963 

7,220 

3,3460 

0,055 

0,1130 

0,0871 

7,158 

3,430 

0,049 

0,112 

14-  15 

0,2253 

1,5820 

0,6790 

0,0060 

0,0050 

0,2037 

1,580 

0,719 

0,006 

0,006 

15-  18 

0,1093 

6,0180 

1,7320 

0,0400 

0,0800 

0,0908 

6,001 

1,794 

0,036 

0,079 

15-  23 

0,1019 

5,0130 

2,9460 

0,0320 

0,0630 

0,0922 

4,957 

2,996 

0,029 

0,063 

16-17 

0,0840 

3,666 

1,4340 

0,0120 

0,0270 

0,0759 

3,609 

1,518 

0,011 

0,027 

18-  19 

0,0651 

2,7780 

0,7520 

0,0050 

0,0100 

0,0589 

2,765 

0,815 

0,005 

0,010 

19-20 

0,0346 

-6,7270 

-2,6580 

0,0170 

0,0340 

0,0313 

-6,740 

-2,595 

0,016 

0,034 

21-22 

0,0118 

-1,871 

-1,593 

0,0010 

0,0001 

0,0106 

-1,899 

-1,599 

0,001 

0,0001 

22-  24 

0,1172 

5,662 

2,787 

0,0440 

0,0670 

0,1060 

5,620 

2,793 

0,039 

0,066 

23-  24 

0,1345 

1,781 

1,282 

0,0060 

0,0120 

0,1217 

1,728 

1,333 

0,005 

0,012 

24  -25 

0,1921 

-1,307 

1,590 

0,0080 

0,0130 

-1,3970 

1,648 

-1,321 

0,008 

0,015 

25  -26 

0,2593 

3,545 

2,367 

0,0450 

0,0670 

0,2345 

3,541 

2,366 

0,041 

0,066 

25-  27 

0,1114 

-4,860 

-0,790 

0,0260 

0,0490 

0,1007 

-4,945 

-0,733 

0,024 

0,050 

27-  28 

0,0000 

-18,173 

-4,169 

2,842  e-14 

1,3090 

0,2952 

-18,227 

-4,110 

-7,270e-ll 

1,312 

27  -29 

0,2240 

6,192 

1,668 

0,0880 

0,1620 

0,000 

6,179 

1,666 

0,079 

0,161 

27-  30 

0,3264 

7,094 

1,662 

0,1650 

0,3040 

0,2026 

7,079 

1,660 

0,148 

0,302 

29-  30 

0,2445 

3,704 

0,606 

0,0340 

0,0630 

0,2211 

-3,669 

-0,542 

0,031 

0,063 

Total  losses 

17,985 

22,407 

Total  losses 

16,079 

21,641 

3.3.  Extreme  learning  machine  algorithm 
3.3.1.  Training-validation  and  test 

The  ELM  training  -validation  and  test  proeess  is  exposed  in  Figure  6,  the  performanee  in  Figure  7,  and  the 
regression  in  Figure  8  respeetively. 
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Figure  6.  Training  and  Testing  Figure  7.  Performance 

Time 

Table  4  illustrates  EFM  training  phase  parameters. 


TaijS 

Figure  8.  Regression 


Table  4.  EFM  training  phase  parameters 


Extreme  learning  maehine 

Values 

Training  time(s) 

0,0469 

Testing  time(s) 

0,0000 

Training  performanee 

3,666e-25 

Testing  performanee 

8,248e-27 

EFM  model  is  built  and  ready  to  be  used  in  real  time. 

Table  5  shows  the  line  flow  results  provided  by  the  EFM  trained  model  for  the  two  resistance  variation: 

Rt=3oc°  and  Rt=o5C° 


Table  5.  Fine  flow  based  on  Extreme  learning  machine  algorithm  in  real  time  (T=3QC°  and  T=5C°) 


i  to  j 

R(pu) 

30C° 

P 

(MW) 

Q 

(MVar) 

Loss 

(MW) 

Loss 

(Mvar) 

R(pu) 

05C° 

P 

(MW) 

Q 

(MVar) 

Loss 

(MW) 

Loss 

(Mvar) 

1-2 

0,0195 

178,010 

-23,312 

5,593 

10,5910 

0,0177 

176,840 

-17,501 

4,960 

10,273 

1-3 

0,0460 

83,370 

4,788 

2,871 

7,1190 

0,0416 

82,639 

6,479 

2,561 

6,961 

2-4 

0,0581 

45,709 

2,525 

1,127 

-0,5200 

0,0525 

45,719 

3,432 

1,024 

-0,507 

2-5 

0,0481 

83,105 

1,320 

3,061 

8,2080 

0,0435 

82,543 

3,232 

2,738 

8,069 

2-6 

0,0592 

61,909 

-1,262 

2,087 

2,2630 

0,0535 

61,918 

0,268 

1,890 

2,270 

3-4 

0,0134 

78,099 

-3,531 

0,788 

1,3500 

0,0121 

77,678 

-1,682 

0,703 

1,320 

4-6 

0,0121 

70,138 

-18,018 

0,618 

1,1880 

0,0109 

70,080 

-15,564 

0,549 

1,147 

4-  12 

0,0000 

44,155 

14,581 

-l,7741e-12 

4,6890 

0,000 

43,989 

14,901 

l,885e-13 

4,669 

5-7 

0,0468 

-14,155 

10,72 

0,156 

-1,6820 

0,0424 

-14,395 

9,613 

0,133 

-1,704 

6-7 

0,0272 

37,485 

-2,102 

0,374 

-0,6000 

0,0246 

37,669 

-1,010 

0,341 

-0,594 

6-8 

0,0122 

29,531 

-4,254 

0,106 

-0,5560 

0,0110 

29,520 

-1,763 

0,094 

-0,563 

6-9 

0,0000 

27,693 

-7,341 

-2,527e-13 

1,5950 

0,0000 

27,690 

-7,249 

-l,175e-13 

1,590 

6-10 

0,0000 

15,822 

0,651 

-4,932e-14 

1,2790 

0,0000 

15,823 

0,658 

-l,779e-14 

1,277 

6-28 

0,0172 

18,809 

-9,686 

0,061 

-13,0820 

0,0155 

18,857 

-9,348 

0,055 

-13,102 

8-  28 

0,0648 

-0,575 

-2,286 

-3,533e-14 

-4,3680 

0,0586 

-0,574 

-2,702 

2,68e-12 

-4,371 

9-  10 

0,0000 

27,693 

6,772 

-2,629e-13 

0,8100 

0,0000 

27,690 

6,617 

-l,848e-14 

0,8060 

9-  11 

0,0000 

-3,388e-13 

-15,708 

-3,024e-14 

0,4650 

0,0000 

2,698e-13 

-15,456 

-2,733e-13 

0,449 

10-17 

0,0330 

5,3610 

4,431 

0,015 

0,0370 

0,0862 

5,415 

4,347 

0,013 

0,037 

10-  20 

0,0954 

9,0270 

3,5720 

0,0820 

0,1810 

0,0298 

9,029 

3,509 

0,074 

0,180 

10-21 

0,0354 

15,7410 

9,8430 

0,1120 

0,2370 

0,0320 

15,701 

9,837 

0,101 

0,236 

10-22 

0,0741 

7,5870 

4,4880 

0,0530 

0,1070 

0,0670 

7,568 

4,499 

0,048 

0,106 

12-  13 

0,0000 

-2,176e-13 

-10,3200 

-l,758e-13 

0,1330 

0,0000 

1.435e-13 

-10.174 

1.691e-13 

0,130 

12-14 

0,1255 

7,8580 

2,433 

0,0760 

0,1550 

0,1135 

7,849 

2,474 

0,069 

0,155 

12-15 

0,0674 

17,8760 

6,9330 

0,2220 

0,4290 

0,0610 

17,783 

7,002 

0,199 

0,426 

12-16 

0,0963 

7,2200 

3,3460 

0,0550 

0,1130 

0,0871 

7,158 

3,430 

0,049 

0,112 

14-  15 

0,2253 

1,5820 

0,6790 

0,0060 

0,0050 

0,2037 

1,580 

0,719 

0,006 

0,006 

15-  18 

0,1093 

6,0180 

1,7320 

0,0400 

0,0800 

0,0908 

6,001 

1,794 

0,036 

0,079 

15-  23 

0,1019 

5,0130 

2,9460 

0,0320 

0,0630 

0,0922 

4,957 

2,996 

0,029 

0,063 

16-17 

0,0840 

3,6660 

1,4340 

0,0120 

0,0270 

0,0759 

3,609 

1,518 

0,011 

0,027 

18-  19 

0,0651 

2,7780 

0,7520 

0,0050 

0,0100 

0,0589 

2,765 

0,815 

0,005 

0,010 

19-20 

0,0346 

-6,7270 

-2,6580 

0,0170 

0,0340 

0,0313 

-6,740 

-2,595 

0,016 

0,034 
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i  to  j 

R(pu) 

30C° 

P 

(MW) 

Q 

(MVar) 

Loss 

(MW) 

Loss 

(Mvar) 

R(pu) 

05C° 

P 

(MW) 

Q 

(MVar) 

Loss 

(MW) 

Loss 

(Mvar) 

21-22 

0,0118 

-1,8710 

-1,5930 

0,0010 

0,0001 

0,0106 

-1,899 

-1,599 

0,001 

0,0001 

22-  24 

0,1172 

5,6620 

2,7870 

0,0440 

0,0670 

0,1060 

5,620 

2,793 

0,039 

0,066 

23-24 

0,1345 

1,7810 

1,2820 

0,0060 

0,0120 

0,1217 

1,728 

1,333 

0,005 

0,012 

24  -25 

0,1921 

-1,3070 

1,5900 

0,0060 

0,0130 

-1,3970 

1,648 

-1,321 

0,008 

0,015 

25  -26 

0,2593 

3,5450 

2,3670 

0,0400 

0,0670 

0,2345 

3,541 

2,366 

0,041 

0,066 

25-27 

0,1114 

-4,8600 

-0,7900 

0,0320 

0,0490 

0,1007 

-4,945 

-0,733 

0,024 

0,050 

27-  28 

0,0000 

-18,173 

-4,1690 

0,0120 

1,3090 

0,2952 

-18,227 

-4,110 

-5,216e-14 

1,312 

27  -29 

0,2240 

6,1920 

1,6680 

0,0050 

0,1620 

0,0000 

6,179 

1,666 

0,079 

0,161 

27-  30 

0,3264 

7,0940 

1,6620 

0,0170 

0,3040 

0,2026 

7,079 

1,660 

0,148 

0,302 

29-  30 

0,2445 

3,7040 

0,6060 

0,0340 

0,0630 

0,2211 

-3,669 

-0,542 

0,031 

0,063 

Total  losses 

17,9850 

22,4070 

Total  losses 

16,079 

21,641 

3.4.  Discussion  of  results 

3.4.1.  Neural  networks 

This  case  study  demonstrates  the  applieation  and  praetieal  benefits  of  the  neural  network  with 
baekpropagation,  as  it  is  introdueed  here  to  solve  a  nonlinear  problem  whieh  is  the  variation  of  the  eonduetor 
resistanee  due  to  the  temperature  fluetuations  and  its  influenee  on  the  power  flow  and  losses  ehanges; 
attempting  to  optimize  the  time  of  simulation. 

The  Newton  raphson  method  generated  the  dataset  to  be  used  for  the  training  patterns  proeess.  The 
input  matrix  is  eomposed  of;  line  resistanee  (R),  load  demand  (Pd),  (QD)(real  and  reaetive  power),  bus 
voltage  magnitude  (Vi).  The  neural  network  model  developed  for  IEEE  30  bus  test  was  trained  to  provide 
power  system  values  whieh  are  power  generation  (real  and  reaetive  power)(Pg),  (Qg)  and  system  losses  (Pij), 
(Qij)(real  and  reaetive). 

The  training  phase  took  15s  to  meet  the  eonvergenee  (Figure  3),  the  performanee  reaehed  the  best  value 
3.15e-19  at  15s  (Figure  4).  The  network  parameters  are  adjusted  (bias  and  weights)  and  it  is  ready  to  be 
utilized  in  real  time  (for  the  two  different  values  of  temperature:  T1=30C°  and  T2=05C°). 

The  eonduetor's  resistanee  inereases  as  its  temperature  inereases  and  so  losses,  due  to;  either  the 
eurrent  load  or  the  weather  eonditions.  Results  shown  in  (Table  3)  are  aequired  in  milliseeonds  (neural 
network  results).  The  objeetive  of  work  is  to  minimize  time  by  utilizing  the  neural  network  for  sueh  ehanging 
of  temperature.  Also  to  assess  the  generated  powers  and  losses  during  the  24  hours  per  day  and  through  the 
entire  year  within  milliseeonds. 

3.4.2.  Extreme  learning  machine 

Most  of  our  effort  has  foeused  on  Extreme  learning  maehine  algorithm  in  this  paper.  Even  for  EEM 
the  simulation  proeess  eontaining  two  steps;  the  training-validation  and  the  test  phase.  The  dataset  is 
randomly  partitioned  into  two  non-overlapped  sets,  to  form  the  input  training  -validation  matrix  and  the  test 
one.  The  input  matrix  is  eomposed  of  line  resistanee  R,  load  demand  (real  and  reaetive  power)  (Pd),  (Qd),  bus 
voltage  magnitude  (Vi).  The  target  matrix  is  eontaining  power  generation  (real  and  reaetive  power))  (Pg),  (Qg) 
and  system  losses  (real  and  reaetive)  (PQ,  (Qij).  The  algorithm  seems  to  be  very  fast,  it  aehieves  the  highest 
aeeuraey  10'^^  (Figure  13)  and  fastest  training  time:  0.0469(S)  and  also  testing  time  0(S)  (figure  12). 

In  order  to  further  verify  the  validity  of  the  improved  built  model,  and  the  generalization;  two  extreme  values 
of  the  resistanee  30C°  and  05C  are  introdueed,  results  are  displayed  in  Table  4. 

Aeeording  to  EEM  based  power  flow,  we  estimated  to  mateh  results  offered  by  the  neural  network 
with  further  optimization  in  time  and  more  exaetitude  (Table  5).  It  is  evident  that  the  regression  aeeuraey  of 
EEM  algorithm  is  superior  to  the  NN. 

Losses  at: 

T=  25C°  are  17,599  MW,  22,244Mvar. 

T=  30C°  are  17,985MW,  22,407  Mvar. 

T=  05C°  are  16,79MW,  21,641Mvar. 

Mathematieally  losses  are  eovered  by  the  plant  eonneeted  to  the  slaek  bus  (slaek  number  one  1),  but 
in  the  real  intereonneeted  power  system;  the  slaek  bus  generator  is  not  always  a  single  equivalent  maehine 
loeated  in  a  eontrol  area,  losses  ean  be  interpreted  as  exehanges  and  power  transaetions  aeross  multi-area 
systems.  This  may  impose  an  additional  operating  eost.  The  priee  inereases  as  well  as  the  slaek  bus 
produetion  inereases. 

ELM  makes  an  improvement  in  the  eleetrieal  system  response;  the  grid  parameters  ean  be  prevented 
instantly  at  any  moment,  and  the  plants  generation  ean  be  reseheduled  to  meet  load  demand  (in  addition  to 
losses),  to  maintain  power  system  seeurity,  to  avoid  eontingeneies  and  load  shedding  (voltage  drop. 
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frequency  fluctuation),  to  prevent  system  violations  and  overloads,  to  predict  the  effect  of  outages  and  to 
minimize  the  fuel  consumption. 

Additionally,  ELM  could  likewise  be  extended  to  make  use  of  a  large  scale  database,  (Within  [05  30]  C°), 
cause  it  is  already  trained,  and  also  for  more  advanced  complications  in  power  system. 

Table  6  lists  the  mean  comparison  results  of  ANN  and  ELM.  Extreme  Learning  Machines  (ELMs) 
are  very  controversial  and  very  fast  machine  learning  models  that  perform  very  well.  ELM  results  are  better 
than  ANN,  due  to  the  ability  of  the  excellent  generalization  performance  and  the  learning  speed  which  can  be 
thousands  of  times  faster  than  traditional  neural  networks. 


Table  6. Comparison 


Neural  network 

Extreme  learning  machine 

Training  time(s) 

15 

0,0469 

Testing  time(s) 

/ 

0,0000 

Training  performance 

3,15e-19 

3,666e-25 

Testing  performance 

/ 

8,248e-27 

4.  CONCLUSION 

Power  system  security  requires  system  monitoring  in  real  time,  where  the  database  is  continuously 
updated.  This  paper  investigated  the  impact  of  the  temperature  variation  on  overhead  power  lines.  This 
influence  on  the  conductor's  resistance  cannot  be  neglected.  We  noticed  different  variations  on  the  behavior 
of  the  electrical  grid,  notably:  the  system  losses,  the  generated  powers,  and  the  line  flow. 

We  introducedExtreme  learning  algorithm  to  solve  this  problem.  The  basic  ANN  modelization  (feed 
forward  back  propagation)  seems  to  be  appropriate  for  the  application  IEEE  30  Bus  system;  but  in  terms  of 
accuracy  and  exactitude  of  results  offered  by  ELM  (the  timing  minimization,  performance),  we  can  say 
economic  objectives  are  realized  (the  optimization  in  time,  fuel  consumption  and  the  reduction  of  the 
operation  cost);  ELM  stands  to  be  the  most  gifted  version  of  the  artificial  intelligence  techniques 
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